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In 2013, a novel betacoronavirus was identified in fecal 
samples from dromedaries in Dubai, United Arab Emirates. 
Antibodies against the recombinant nucleocapsid protein of 
the virus, which we named dromedary camel coronavirus 
(DcCoV) UAE-HKU23, were detected in 52% of 59 drom- 
edary serum samples tested. In an analysis of 3 complete 
DcCoV UAE-HKU23 genomes, we identified the virus as 
a betacoronavirus in lineage A1. The DcCoV UAE-HKU23 
genome has G+C contents; a general preference for G/C in 
the third position of codons; a cleavage site for spike pro- 
tein; and a membrane protein of similar length to that of 
other betacoronavirus A1 members, to which DcCoV UAE- 
HKU23 is phylogenetically closely related. Along with this 
coronavirus, viruses of at least 8 other families have been 
found to infect camels. Because camels have a close as- 
sociation with humans, continuous surveillance should be 
conducted to understand the potential for virus emergence 
in camels and for virus transmission to humans. 


Ihe 2003 epidemic of severe acute respiratory syn- 

drome (SARS) boosted interest in the discovery of 
new coronaviruses (CoVs) (/—3). In 2004, a novel human 
CoV (HCoV), named HCoV-NL63, was reported (4), and 
the discovery of another novel HCoV, HCoV-HKU1, was 
described and further characterized in 2005 (5,6) and 2006 
(7). SARS-CoV-like viruses have also been reported in 
Chinese horseshoe bats in Hong Kong, China, and other 
horseshoe bats in China (8,9). The discovery in Chinese 
horseshoe bats in Yunnan, China, of a new SARS-CoV— 
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like virus that uses ACE2 as receptor has furthered inter- 
est in discovering animal origins of human infections (/0). 
We have discovered 20 other animal CoVs that include 2 
novel betacoronavitus lineages and a novel genus, Delta- 
coronavirus (11—20). From our studies, bats and birds were 
shown to be the gene sources for fueling the evolution and 
dissemination of alphacoronaviruses and betacoronavi- 
ruses and of gammacoronaviruses and deltacoronaviruses, 
respectively (/8). 

In 2012, a novel CoV, Middle East respiratory syn- 
drome CoV (MERS-CoV) emerged as a cause of severe 
respiratory infections associated with high rates of death 
among humans; the virus is closely related to tylonycteris 
bat CoV HKU4 and pipistrellus bat CoV HKUS (Pi-Bat 
CoV HKUS) (2/—23). It has also been shown that drome- 
daries in the Middle East possess MERS-CoV neutralizing 
antibodies (24). To further knowledge of the evolution and 
dissemination of CoVs, we conducted a molecular epide- 
miology study of fecal samples obtained from dromedaries 
in Dubai, United Arab Emirates. 


Materials and Methods 


Samples 

Dromedary fecal and serum samples used in the study 
were leftover specimens that had been submitted for patho- 
gen screening (feces) or preventive health screening (se- 
rum) to Central Veterinary Research Laboratory (Dubai, 
United Arab Emirates) during January—July 2013. The fe- 
cal and serum samples were not obtained from the same 
animals. None of the dromedaries tested were known to 
have had contact with bats or horses. 

We tested a total of 293 fecal samples: 232 from teen- 
age and adult dromedaries (Camelus dromedarius) (>1 
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year of age) and 61 from dromedary calves (<1 year of 
age). Among the 293 samples, 6 were collected in January, 
209 in February, 5 in March, 39 in April, 16 in May, 7 in 
June, and 11 in July 2013. 

We tested a total of 59 serum samples: 55 from teen- 
age and 4 from adult dromedaries. The serum samples were 
collected from female dairy farm or racing dromedaries. 
The dairy dromedaries were purchased from various coun- 
tries (e.g., Saudi Arabia, Oman, Sudan, and Pakistan); the 
number of dairy dromedaries from each country and their 
lengths of stay in Dubai were not known. The racing drom- 
edaries were from Dubai Emirate. 


RNA Extraction, Reverse Transcription PCR, 
and DNA Sequencing 

Viral RNA extraction was conducted as described 
(20,23). Initial CoV screening was performed by ampli- 
fying a 440-bp fragment of CoV RNA-dependent RNA 
polymerase (RdRp) gene by using conserved primers 
(5'-GGTTGGGACTATCCTAAGTGTGA-3’ and 5’-AC- 
CATCATCNGANARDATCATNA-3’ and degenerative 
bases N, R, and D, representing A/C/T/G, A/G, and A/G/T, 
respectively. After the novel CoV was detected in samples, 
subsequent screening was performed by using specific prim- 
ers 5'-ACTATGACTGGCAGAATGTT-3’ and 5'-TAATA- 
AGGCGACGTAACATA-3’. To amplify a 126-bp fragment 
of RdRp gene, reverse transcription PCR (RT-PCR) and 
DNA sequencing were performed as described (20,23). 


Virus Culture 

The fecal samples from 3 dromedaries tested positive 
for CoV. These samples were cultured in HRT-18G, Vero 
E6, Caco-2, and LLC-MK2 cell lines. 


Complete Genome Sequencing and Analysis 

Three complete genomes of DcCoV UAE-HKU23 
(265F, 362F, and 368F) were amplified and sequenced 
as described (5,/2). RNA extracted from fecal specimens 
was used as template, and a database of CoV genes and 
genomes (CoVDB, http://covdb.microbiology.hku.hk) (25) 
was used for sequence retrieval. Sequences were assem- 
bled and edited to produce final sequences. 

We used EMBOSS Needle (http://www.ebi.ac.uk/ 
Tools/psa/emboss_needle/) to compare the nucleotide se- 
quences of the genomes and the deduced amino acid se- 
quences with those for other CoVs. A neighbor-joining 
phylogenetic tree with 1,000 bootstraps was constructed by 
using the Jones-Taylor-Thornton substitution model; gam- 
ma distribution among sites was conducted in MEGAS (26). 


Real-time Quantitative RT-PCR 
According to our protocol (19,20), all samples 
positive for DcCoV UAE-HKU23 by RT-PCR were 
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subjected to quantitative RT-PCR (qRT-PCR).The as- 
say was performed, as described (27), using a real- 
time l-step qRT-PCR with DcCoV UAE-HKU23 
primers 5'-ATAGCGGCTACACGTGGTGTT-3' and 
5'-TCCCAGCCGCCATAAAACT-3’ and probe 5’- 
(FAM) CTGTTGTTATAGGCACCACT (BHQ1)-3". To 
generate calibration curves, we prepared a series of 6 log,, 
dilutions equivalent to 10'-10° copies per reaction mix- 
ture and ran them in parallel with the test samples. 


Cloning and Purification 

The nucleocapsid proteins of SARS-CoV (betacorona- 
virus lineage B), Pi-Bat CoV HKUS (betacoronavirus lin- 
eage C), and rousettus bat coronavirus (Ro-BatCov HKU9; 
betacoronavirus lineage D) were obtained as described 
(2,16). To produce a plasmid for DceCoV UAE-HKU23 
nucleocapsid protein purification, primers 5’-CATGC- 
CATGGGCATGTCTTTTACTCCTGGTAAGC-3’ and 5’- 
CCGCTCGAGTATTTCTGAGGTGTTTTCAG-3’ were 
used to amplify the gene encoding the nucleocapsid protein 
of DcCoV UAE-HKU23. The sequence encoding aa 1-672 
of the nucleocapsid protein was amplified and cloned into 
the NcoI and Xhol sites of expression vector pET-28b(+) 
(Merck KGaA, Darmstadt, Germany). Recombinant nucleo- 
capsid protein was expressed and purified by using the Ni?*- 
loaded HiTrap Chelating System (GE Healthcare Life Sci- 
ences, Little Chalfont, UK) according to the manufacturer’s 
instructions. The nucleocapsid protein of MERS-CoV was 
cloned and purified by using the method described above 
with primers 5'-GGAATTCCATATGATGGCATCCCCT- 
GCTGCACCTC-3' and 5'-ATAAGAATGCGGCCGCAT- 
CAGTGTTAACATCAATCATT-3’. 


Western Blot Analysis 

Western blot analysis was performed as described 
(5) with 1.5 ug of purified (His),-tagged recombinant nu- 
cleocapsid protein of DcCoV UAE-HKU23 and 1:2,000, 
1:4,000, and 1:8,000 dilutions of dromedary serum samples. 
Antigen—antibody interaction was detected by using 1:4,000 
diluted horseradish peroxidase—conjugated Goat Anti-Lla- 
ma IgG (Life Technologies, Carlsbad, CA, USA) and the 
ECL fluorescence system (GE Healthcare Life Sciences). 

To determine the possibility of cross-reactivity be- 
tween antibodies against the nucleocapsid protein of 
DcCoV UAE-HKU23 and that of other betacoronavirus 
lineages, we tested 3 serum samples that were positive 
for antibody against the nucleocapsid protein of Dc- 
CoV UAE-HKU23. We used 1.5 wg of purified (His),- 
tagged recombinant nucleocapsid protein of 3 betacoro- 
naviruses (SARS-CoV [lineage B], Pi-BatCoV HKUS5 
[lineage C], and Ro-BatCoV HKU9 [lineage D]) and 
1:2,000 dilutions of serum samples. To determine the 
presence of antibodies against the nucleocapsid protein 


561 


RESEARCH 


of MERS-CoV, we tested the serum samples by using 1.5 
ug of purified (His),.-tagged recombinant nucleocapsid 
protein of MERS-CoV and 1:2,000 dilutions of serum 
samples. Antigen—antibody interaction was detected as 
described above. 


Indirect Immunofluorescence 

Anti-MERS-CoV antibody detection by indirect im- 
munofluorescence was performed as described (28) with 
minor modifications. In brief, Vero cells infected with 
MERS-CoV were prepared as described (28). Camel serum 
samples were screened at a dilution of 1:160 on infected 
and noninfected control cells. Antigen—antibody interac- 
tion was detected by using fluorescein isothiocyanate—con- 
jugated Goat Anti-Llama IgG (Life Technologies). Serum 
samples positive at a screening dilution of 1:160 were fur- 
ther titrated with serial 2-fold dilutions. The indirect im- 
munofluorescence antibody titer was the highest dilution 
giving a positive result. 


Neutralization Antibody Test 

The neutralization antibody test was performed as 
described (28). In brief, starting with a serum dilution of 
1:10, we prepared serial 2-fold dilutions of serum in 96- 
well microtiter plates. For each serum dilution, 0.05 mL of 
the dilution was mixed with 0.05 mL of 200 MERS-CoV 
50% tissue culture infectious doses and incubated at 37°C 
for 1.5 h in a CO, incubator. Then 0.1 mL of virus—serum 
mixture was inoculated in duplicate wells of 96-well mi- 
crotiter plates with preformed monolayers of Vero cells 
and further incubated at 37°C for 3-4 days. Cytopathic ef- 
fects were observed by using an inverted microscope on 
days 3 and 4 after inoculation. The neutralizing antibody 
titer was determined as the highest dilution of serum that 
completely suppressed the cytopathic effects in at least half 
of the infected wells. 


Estimation of Substitution Rates and 
Divergence Dates 

The number of synonymous substitutions per synony- 
mous site, Ks, and the number of nonsynonymous substitu- 
tions per nonsynonymous site, Ka, for each coding region 
between each pair of strains were calculated by using the 
Nei-Gojobori method (Jukes-Cantor) in MEGAS (26). Di- 
vergence time was calculated on the basis of RdRp gene 
sequence data by using a Bayesian Markov Chain Monte 
Carlo approach as implemented in BEAST version 1.8.0 
(http://beast.bio.ed.ac.uk), as described (/5,/9,29,30). 
Bayesian skyline under a relaxed-clock model with an un- 
correlated exponential distribution was adopted for making 
inferences because Bayes factor analysis for the RdRp gene 
indicated that this model fitted the data better than other 
models tested. 
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Identification of CoV in Dromedaries 

Of the 293 fecal samples tested, 14 (4.8%) were RT- 
PCR positive for the CoV RdRp gene; 1 (0.4%) of the 
samples was from an adult dromedary and 13 (21.3%) 
were from calves (Table 1). Of the 14 positive samples, 11 
were collected in April and 3 were collected in May. Ten 
of the 14 samples were submitted to the Central Veterinary 
Research Laboratory for routine checking, and the other 4 
were collected because the dromedaries had diarrhea. Se- 
quencing results indicated a 126-nt sequence identical to 
that of equine CoV, a betacoronavirus | species in lineage 
A (betacoronavirus A 1). 


Virus Culture and Virus Load 

Attempts to stably passage DcCoV UAE-HKU23 in 
cell cultures were unsuccessful; no cytopathic effect or vi- 
ral replication was detected. Real-time qRT-PCR showed 
that the amounts of DcCoV UAE-HKU23 RNA ranged 
from 5.7 x 10* copies/mL to 9.8 x 10’ copies/mL (medi- 
an 8.4 x 10°) in the 14 fecal samples positive for DcCoV 
UAE-HKU23 (Table 1). 


Genome Organization and Coding Potential 

The 3 complete genomes of DcCoV UAE-HKU23 
(GenBank accession nos. KF906249-KF906251) were 
31,036 bases and had a G+C content of 37% (Table 2). 
The genome organization is similar to that of other beta- 
coronavirus lineage A CoVs (Figure 1). Additional open- 
reading frames (ORFs) coding for nonstructural proteins 
(NSPs) NS2 and NSS are found. DeCoV UAE-HKU23 
and other CoVs in betacoronavirus lineage A possess the 
same putative transcription regulatory sequence motif, 
5'-UCUAAAC-3’, at the 3’ end of the leader sequence 
and preceding most ORFs (Table 3; online Technical 


Table 1. Epidemiologic data for dromedaries in the Middle East 
that were positive for a novel betacoronavirus, DcCoV UAE- 
HKU23, 2013* 


Age category Reason for sample testing Virus load 
Calf Routine check 3.3 x 10° 
Calf Routine check 1.1 x 10° 
Calf Routine check 1.5 x 10° 
Calf Routine check 3.4 x 10° 
Calf Routine check 3.2 x 10° 
Calf Routine check 6.5 x 10° 
Calf Routine check 9.5 x 10° 
Calf Routine check 5.7 x 10° 
Calf Routine check 2.1 x 10° 
Calf Routine check 4.4 x 10° 
Adult Diarrhea 7.3 x 10° 
Calf Diarrhea 9.8 x 10” 
Calf Diarrhea 4.4 x 10” 
Calf Diarrhea 1.0 x 10’ 


*DcCoV, dromedary camel coronavirus; calf, <1 year of age. 
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Table 2. Comparison of representative coronaviruses with a novel betacoronavirus, DcCoV UAE-HKU23, discovered in dromedaries in 


the Middle East, 2013* 


Genome features 


% Pairwise amino acid identity with DcCoV UAE-HKU23 (265F) 


Coronavirus (genome) Size, bases | G+C content 3CL"° RdRp Helicase S protein N protein 
Alphacoronavirus 
HCoV-229E 27,317 0.38 44.4 55.8 57.9 26.9 25.1 
HCoV-NL63 27,553 0.34 42.8 56.0 57.5 26.4 24.1 
Betacoronavirus 
Lineage A 
BCoV 31,028 0.37 99.7 99.9 99.5 94.1 98.4 
SACoV 30,995 0.37 99.7 99.8 99.3 94.1 98.9 
DcCoV UAE-HKU23 (362F) 31,036 0.37 99.7 100 100 99.8 100 
DcCoV UAE-HKU23 (368F) 31,036 0.37 99.7 100 100 99.8 100 
Lineage B 
SARS-CoV 29,751 0.41 48.4 66.3 68.0 30.4 34.4 
Lineage C 
MERS-CoV 30,119 0.41 53.3 68.2 67.9 31.8 34.7 
Lineage D 
Ro-BatCoV HKU9 29,114 0.41 48.4 66.5 67.5 29.2 32.5 
Gammacoronavirus 
Infectious bronchitis virus 27,608 0.38 42.1 60.6 60.2 24.8 28.0 
BdCoV HKU22 31,750 0.39 45.3 60.0 58.2 25.7 28.9 
Deltacoronavirus 
BuCoV HKU11 26,476 0.39 38.5 51.2 47.5 26.4 23.8 
PorCoV HKU15 25,421 0.43 38.8 51.8 48.6 26.0 23.1 


*The representative coronaviruses were those with complete genomes sequences available. DcCoV, dromedary camel coronavirus; 3CL"°, chymotrypsin- 
like protease; RdRp, RNA-dependent RNA polymerase; S protein, spike protein; N protein, nucleocapsid protein; HCoV, human CoV; BCoV, bovine CoV; 
SACoV, sable antelope CoV; SARS-CoV, severe acute respiratory syndrome CoV; MERS-CoV, Middle East respiratory syndrome CoV; Ro-BatCoV, 
rousettus bat CoV; BdCoV, bottlenose dolphin CoV; BuCoV, bulbul CoV; PorCoV, porcine CoV. 


Appendix Table 1, wwwnc.cdc.gov/EID/article/20/4/13- 
1769-Techapp1.pdf) (19,30). 

The characteristics of putative NSPs in ORFlab of 
DcCoV UAE-HKU23 are shown in online Technical Ap- 
pendix Table 2. The ORF lab polyprotein shared 70.7%-— 
99.3% aa identity with polyproteins of other betacoronavi- 
tus lineage A CoVs. The predicted putative cleavage sites 
were conserved between DcCoV UAE-HKU23 and other 
members of betacoronavirus Al of the Betacoronavirus ge- 
nus. The lengths of NSPs 1-3, 13, and 15 in DcCoV UAE- 
HKU23 differed from those in equine CoV, porcine hem- 
agglutinating encephalomyelitis virus, and/or HCoV-OC43 
as a result of deletions/insertions. 

The amino acid sequence of the predicted spike pro- 
tein of DcCoV UAE-HKU23 is most similar to that of 
bovine coronavirus (BCoV) and sable antelope CoV, 
with which DcCoV UAE-HKU23 has 94.1% similarity 
(Table 2). A comparison of the amino acid sequences of 
DcCoV UAE-HKU23 spike protein and BCoV spike pro- 
tein showed 81 aa polymorphisms, of which 24 were seen 
within the region previously identified as hypervariable 
among the spike protein of other betacoronavirus lineage 
A CoVs (3/) (Figure 2); this finding suggests that this 
region in DcCoV UAE-HKU23 is also subject to strong 
immune selection. BCoV has been found to use N-acetyl- 
9-O acetyl neuramic acid as a receptor for initiation of in- 
fection (32). Among the 5 aa that may affect S1-mediated 
receptor binding in BCoV (3/), 2 aa (threonine at posi- 
tion 11 and glutamine at position 179) were conserved in 
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DcCoV UAE-HKU23 (Figure 2). However, at positions 
115, 118, and 173, aspartic acid, methionine, and as- 
paragine observed in BCoV and were replaced by serine, 
threonine, and histidine, respectively, in DcCoV UAE- 
HKU23. A recent report identified 4 aa acids that were 
critical sugar-binding residues in the spike protein of 
BCoV (tyrosine, glutamic acid, tryptophan, and histidine 
at positions 162, 182, 184, and 185, respectively) (32); all 
4 aa were also present in spike protein of DcCoV UAE- 
HKU23 (Figure 2). Another study identified 7 aa substitu- 
tions in the spike protein of BCoV that differed between 
virulent and avirulent, cell culture-adapted strains (3/); 5 
of the 7 aa from virulent strains were also conserved in 
DcCoV UAE-HKU23, and amino acid substitutions were 
observed in the other 2 aa (threonine—valine at position 
40 and aspartic acid—asparagine at position 470). It has 
also been reported that an amino acid change at position 
531 of the spike protein of BCoV discriminated between 
enteric (aspartic acid/asparagine) and respiratory (gly- 
cine) strains (33). At this position, a threonine was con- 
served in all 3 genomes of DcCoV UAE-HKU23. 

NSS of DcCoV UAE-HKU23 shares 83.5%-98.2% 
aa identity with the corresponding NSPs of betacoronavi- 
rus Al members. In murine hepatitis virus, translation of 
the envelope protein is cap-independent, through an in- 
ternal ribosomal entry site (19,30). However, a preceding 
transcription regulatory sequence, 5'’-UCCAAAC-3’, can 
be identified upstream of the envelope protein in DcCoV 
UAE-HKU23, as in other betacoronavirus Al members 
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Figure 1. Genome organizations of a novel betacoronavirus, in boldface, discovered in dromedaries in the Middle East in 2013, and 
representative coronaviruses from each coronavirus genus (labeled on left). Numbers at top represent genome position. A, B, C, and 
D represent betacoronavirus lineages. Papain-like proteases (PL1°, PL2°°, and PL*°), chymotrypsin-like protease (3CL°°), and RNA- 
dependent RNA polymerase (RdRp) are represented by orange boxes. Hemagglutinin esterase (HE), spike (S), envelope (E), membrane 
(M), and nucleocapsid (N) proteins are represented by green boxes. Putative accessory proteins are represented by blue boxes. HCoV, 
human coronavirus; BCoV, bovine CoV; DcCoV, dromedary camel CoV; SARS-CoV, severe acute respiratory syndrome CoV; MERS-CoV, 
Middle East respiratory syndrome CoV; Ro-BatCoV, rousettus bat CoV; IBV, infectious bronchitis virus; BdCoV, bottlenose dolphin CoV; 
BuCoV, bulbul CoV; PorCoV, porcine CoV. 


(Table 3; online Technical Appendix Table 1) (/9,30). 
Downstream to nucleocapsid protein, the 3’-untranslated 
region contains a predicted bulged stem-loop structure 
of 68 nt (nt position 30747-30814) that is conserved in 
betacoronaviruses (34). Overlapping with the bulged 
stem-loop structure by 5 nt, is a conserved pseudoknot 


structure (nt position 30810-30863) that is important for 
CoV replication. 


Phylogenetic Analyses 
Phylogenetic trees constructed by using the amino 
acid sequences of ORF 1b polyprotein, spike protein, and 


Table 3. Putative transcription regulatory sequence of betacoronavirus A1 members* 


Putative transcription regulatory sequence (distance, in bases, to AUG)t 


DcCoV UAE- 
ORF _HKU23 BCoV CRCoV SACoV GiCoV PHEV HCoV-OC43 ECoV 
Jab UCUAAAC UCUAAAC  UCUAAAC  UCUAAAC  UCUAAAC  UCUAAAC  UCUAAAC  UCUAAAC (140) 
(140) (140) (140) (140) (140) (140) (140) 
NS2 UCUAAAC UCUAAAC UCUAAAC UCUAAAC UCUAAAC  UCUAAAC ~=UCUAAAC —sUCUAAAA (12) 
(7) (7) (7) (7) (7) (1) (7) 
HE | ACUAAAC ACUAAAC ACUAAAC ACUAAAC ACUAAAC ACUAAAC = AUUAAAC UCUAAAC (9) 
(9) (9) (9) (9) (9) (9) (9) 
iS UCUAAAC UCUAAAC UCUAAAC UCUAAAC UCUAAAC UCUAAAC  UCUAAAC _—_—UCUAAAC (0) 
(0) (0) (0) (0) (0) (0) (0) 
NS5 GGUGAAC GGUAGAC GGUAGAC GGUAGAC GGUAGAC UUAAGCA UCUAGCA UAUACUUUAUAA 
(51) (50) (51) (50) (50) (32) (20/32) (41) 
E UCCAAAC UCCAAAC UCCAAAC UCCAAAC UCCAAAC UCCAAAC UCCAAAC — UCCAAAC (123) 
(123) (123) (123) (123) (123) (123) (123) 
M UCCAAAC UCCAAAC UCCAAAC UCCAAAC UCCAAAC UCCAAAC  UCCAAAC __UCCAAAC (3) 
(3) (3) (3) (3) (3) (3) (3) 
N UCUAAAC UCUAAAC UCUAAAC  UCUAAAC  ~=UCUAAAC  ~=UCUAAAC. ~—s UCUAAAU UCUAAAC (7) 


(7) (7) (7) (7) 


(7) (7) (7) 


*ORF, open reading frame; DcCoV, dromedary camel coronavirus; BCoV, bovine CoV; CRCoV, canine respiratory CoV; SACoV, sable antelopeCoV; 
GiCoV, giraffe CoV; PHEV, porcine hemagglutinating encephalomyelitis virus; HCoV, human CoV; ECoV, equine CoV; NS, nonstructural protein; HE, 
hemagglutinin esterase; S, spike protein; E, envelope protein; M, membrane protein; N, nucleocapsid protein. 

tBoldface indicates bases and distances to AUG that are different from those for DcCoV UAE-HKU23. 
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Figure 2. Amino acid comparison of the spike protein of a novel betacoronavirus, dromedary camel coronavirus (DcCCoV) UAE-HKU23, 
discovered in dromedaries in the Middle East in 2013, with that of bovine coronavirus (BCoV; GenBank accession no. AF391541). Amino 
acid substitution sites, key amino acids for virulence and receptor binding in BCoV, and cleavage sites are shown. Amino acid positions 
are given with reference to DcCoV UAE-HKU23. Conserved amino acids, compared with those of DcCoV UAE-HKU23 (strain 265F), 
are represented by dots. Amino acids of putative cleavage sites are underlined. Amino acids within the S1 hypervariable region of BCoV 
are marked with open boxes. Amino acid sites central to virulence in BCoV are highlighted in light gray. Amino acid sites shown to affect 
$1-mediated receptor binding in BCoV are highlighted in dark gray. The 4 critical sugar-binding residues are indicated by black dots. The 
amino acid site that discriminated between enteric and respiratory BCoV strains is highlighted in black. 


nucleocapsid protein of DcCoV UAE-HKU23 and other 
CoVs are shown in Figures 3—5. The pairwise amino acid 
identities of chymotrypsin-like protease (3CL"”), RdRp, 
helicase, spike protein, and nucleocapsid protein are shown 
in Table 2. In all 3 phylogenetic trees, DcCoV UAE- 
HKU23 clustered with other members of betacoronavirus 
Al, including BCoV, sable antelope CoV, equine CoV, 
HCoV-O0C43, giraffe CoV, porcine hemagglutinating en- 
cephalomyelitis virus, canine respiratory CoV, and rabbit 
CoV (RbCoV) HKU14 (Figures 3-5). 


Antibody Detection 

Nucleocapsid protein of DcCoV UAE-HKU23 was 
purified. Prominent immunoreactive bands were visible for 
31 (52%) of 59 dromedary serum samples; 25 of the 31 
samples had titers of 2,000, three had titers of 4,000, and 
3 had titers of 8,000 (Figure 6). Serum samples for all 4 
adult dromedaries were positive for DoCoV UAE-HKU23 
antibodies, and 27 (49%) of the 55 samples for teenage 
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dromedaries were positive. Band sizes were ~50 kDa, con- 
sistent with the expected size of 50.4 kDa for the full-length 
(His),-tagged recombinant nucleocapsid protein. Only very 
faint bands were observed when the 3 serum samples posi- 
tive for DcCoV UAE-HKU23 antibodies were incubated 
with nucleocapsid proteins of SARS-CoV, Pi-BatCoV 
HKUS, or Ro-BatCoV HKU9, indicating minimal cross- 
reactivity. This finding concurs with our previous observa- 
tion that minimal cross-reactivity occurs between CoVs in 
different lineages in betacoronavirus (/6). For MERS-CoV 
antibody testing, results were positive for 57 (97%) of the 
59 samples by Western blot analysis, for all 59 samples by 
indirect immunofluorescence, and for 58 (98%) of the 59 
samples by neutralization antibody test (Table 4). 


Estimation of Substitution Rates and 
Divergence Dates 

The Ka, Ks, and Ka/Ks of the various coding re- 
gions in DcCoV UAE-HKU23 are shown in Table 5. 
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Figure 3. Phylogenetic analysis of open reading frame (ORF) 1b polyprotein of dromedary camel coronavirus (DcCoV) UAE-HKU23 from 
dromedaries of the Middle East, 2013. The tree was constructed by the neighbor-joining method, using Jones-Taylor-Thornton substitution 
model with gamma distributed rate variation and bootstrap values calculated from 1,000 trees. Bootstrap values of <70% are not shown. 
A total of 2,725 aa positions in ORF 1b polyprotein were included in the analysis. The tree was rooted to Breda virus (GenBank accession 
no. AY_427798). Betacoronavirus lineages A1 and A-D are indicated on the right. Boldface indicates the 3 strains of DCCoV UAE-HKU23 
characterized in this study. Virus definitions and GenBank accession numbers (in parentheses) follow: Mi-BatCoV 1B, miniopterus bat 
CoV 1B (NC_010436); Mi-BatCoV 1A (NC_010437); Mi-BatCoV HKU8 (NC_010438); Ro-BatCoV HKU10, rousettus bat CoV HKU10 
(JQ989270); Hi-BatCoV HKU10, hipposideros bat CoV HKU10 (JQ989266); HCoV-NL63, human CoV NL63 (NC_005831); HCoV-229E 
(NC_002645); RhBatCoV HKU2, rhinolophus bat CoV HKU2 (EF203064); Sc-BatCoV-512, scotophilus bat CoV 512 (NC_009657); PEDV, 
porcine epidemic diarrhea virus (NC_003436); MCoV, mink CoV (HM245925); CCoV, canine CoV (GQ477367); FIPV, feline infectious 
peritonitis virus (AY994055); PRCV, porcine respiratory CoV; TGEV, transmissible gastroenteritis virus; SACoV, sable antelope CoV; 
GiCoV, giraffe CoV (EF424622); CRCoV, canine respiratory CoV (JX860640); BCoV, bovine CoV (NC_003045); HCoV-OC43, human CoV 
OC43 (NC_005147); PHEV, porcine hemagglutinating encephalomyelitis virus (NC_007732); ECoV, equine CoV (NC_010327); RbCoV 
HKU14, rabbit CoV HKU14 (JN874559); HCoV-HKU1, human CoV HKU1 (NC_006577); MHV, murine hepatitis virus (NC_001846); 
RCoV, rat CoV (NC_012936); Pi-BatCoV HKUS5, pipistrellus bat CoV HKU5 (NC_009020); MERS-CoV, Middle East respiratory syndrome 
CoV (JX869059); Ty-BatCoV HKU4, tylonycteris bat CoV HKU4 (NC_009019); SARSr-CiCoV, SARS-related palm civet CoV (AY304488); 
SARS-CoV, severe acute respiratory syndrome—associated human CoV (NC_004718); SARSrCoV CFB, SARS-related Chinese 
ferret badger CoV (AY545919); SARS-related rhinolophus bat CoV HKU3 (DQ022305); Ro-BatCoV HKU9, rousettus bat CoV HKU9 
(NC_009021); TCoV, turkey CoV (NC_010800); IBV-peafowl, peafowl CoV (AY641576); IBV, infectious bronchitis virus (NC_001451); 
IBV-partridge, partridge CoV (AY646283); BdCoV HKU22, bottlenose dolphin CoV HKU22 (KF793824); BWCoV-SW1, Beluga whale 
CoV SW1 (NC_010646); WiCoV HKU20, wigeon CoV HKU20 (JQ065048); NHCoV HKU19, night-heron CoV HKU19 (JQ065047); 
CMCoV HKU21, common-moorhen CoV HKU21 (JQ065049); BuCoV HKU11, bulbul CoV HKU11 (FJ376619); ThCoV HKU12, thrush 
CoV HKU 12 (FJ376621); WECoV HKU16, white-eye CoV HKU16 (JQ065044); MunCoV HKU13, munia CoV HKU13 (FJ376622); MRCoV 
HKU18, magpie—-robin CoV HKU18 (JQ065046); PorCoV HKU15, porcine CoV HKU15 (JQ065042):; SpCoV HKU17, sparrow CoV HKU17 
(JQ065045). Numbers at nodes represent bootstrap values. Scale bar indicates the estimated number of substitutions per 10 aa. 
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Figure 4. Phylogenetic analyses of spike protein of dromedary camel coronavirus (DcCoV) UAE-HKU23 from dromedaries of the Middle 
East, 2013. The tree was constructed by the neighbor-joining method, using Jones-Taylor-Thornton substitution model with gamma 
distributed rate variation and bootstrap values calculated from 1,000 trees. Bootstrap values of <70% are not shown. A total of 1,366 
aa positions in spike protein were included in the analysis. The tree was rooted to Breda virus (GenBank accession no. AY_427798). 
Betacoronavirus lineages A1 and A-D are indicated on the right. Boldface indicates the 3 strains of DcCoV UAE-HKU23 characterized in 
this study. Virus definitions and GenBank accession numbers (in parentheses) follow: RbCoV HKU14, rabbit CoV HKU14 (JN874559); 
HCoV-0C43, human CoV OC43 (NC_005147); CRCoV, canine respiratory CoV (JX860640); BCoV, bovine CoV (NC_003045); SACoV, 
sable antelope CoV (EF424621); GiCoV, giraffe CoV (EF424622); PHEV, porcine hemagglutinating encephalomyelitis virus (NC_007732); 
ECoV, equine CoV (NC_010327); MHV, murine hepatitis virus (NC_001846); RCovV, rat CoV (NC_012936); HCoV-HKU1, human CoV 
HKU1 (NC_006577); MERS-CoV, Middle East respiratory syndrome CoV (JX869059); Ty-BatCoV HKU4, tylonycteris bat CoV HKU4 
(NC_009019); Pi-BatCoV HKU5, pipistrellus bat CoV HKU5 (NC_009020); SARSrCoV CFB, SARS-related Chinese ferret badger CoV 
(AY545919); SARSr-CiCoV, SARS-related palm civet CoV (AY304488); SARS-CoV, severe acute respiratory syndrome—associated 
human CoV (NC_004718); SARSr-Rh-BatCoV HKU3, SARS-related rhinolophus bat CoV HKU3 (DQ022305); Ro-BatCoV HKU9, 
rousettus bat CoV HKU9 (NC_009021); RhBatCoV HKU2, rhinolophus bat CoV HKU2 (EF203064); BWCoV-SW1, Beluga whale CoV 
SW1 (NC_010646); BdCoV HKU22, bottlenose dolphin CoV HKU22 (KF793824); TCoV, turkey CoV (NC_010800); IBV-partridge, partridge 
CoV (AY646283); IBV, infectious bronchitis virus (NC_001451); IBV-peafowl, peafowl CoV (AY641576); CCoV, canine CoV (GQ477367); 
FIPV, feline infectious peritonitis virus (AY994055); PRCV, porcine respiratory CoV (DQ811787); TGEV, transmissible gastroenteritis virus 
(DQ811789); MCoV, mink CoV (HM245925); Sc-BatCoV-512, scotophilus bat CoV 512 (NC_009657); PEDV, porcine epidemic diarrhea 
virus (NC_003436); Mi-BatCoV 1B, miniopterus bat CoV 1B (NC_010436); Mi-BatCoV 1A, miniopterus bat CoV 1A (NC_010437); Mi- 
BatCoV HKU8, miniopterus bat CoV HKU8 (NC_010438); Ro-BatCoV HKU10, rousettus bat CoV HKU10 (JQ989270); Hi-BatCoV HKU10, 
hipposideros bat CoV HKU10 (JQ989266); HCoV-NL63, human CoV NL63 (NC_005831); HCoV-229E, human CoV 229E (NC_002645); 
WiCoV HKU20, wigeon CoV HKU20 (JQ065048); MRCoV HKU18, magpie—robin CoV HKU18 (JQ065046); SoCoV HKU17, sparrow CoV 
HKU17 (JQ065045); NHCoV HKU19, night-heron CoV HKU19 (JQ065047); ThCoV HKU12, thrush CoV HKU12 (FJ376621); CMCoV 
HKU21, common-moorhen CoV HKU21 (JQ065049); WECoV HKU16, white-eye CoV HKU16 (JQ065044); MunCoV HKU13, munia CoV 
HKU13 (FJ376622); PorCoV HKU15, porcine CoV HKU15 (JQ065042); BuCoV HKU11, bulbul CoV HKU11 (FJ376619). Numbers at 
nodes represent bootstrap values. Scale bar indicates the estimated number of substitutions per 5 aa. 
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Figure 5. Phylogenetic analyses of the nucleocapsid protein of a novel coronavirus (CoV), dromedary camel CoV (DcCoV) UAE-HKU23, 
discovered in dromedaries of the Middle East, 2013. The tree was constructed by the neighbor-joining method, using Jones-Taylor- 
Thornton substitution model with gamma distributed rate variation and bootstrap values calculated from 1,000 trees. Bootstrap values of 
<70% are not shown. A total of 448 aa positions were included in the analysis. The tree was rooted to Breda virus (GenBank accession 
no. AY_427798). Betacoronavirus lineages A1 and A-D are indicated on the right. Boldface indicates the 3 strains of DcCoV UAE- 
HKU23 characterized in this study. Virus definitions and GenBank accession numbers (in parentheses) follow: TGEV, transmissible 
gastroenteritis virus (DQ811789); PRCV, porcine respiratory CoV (DQ811787); CCoV, canine CoV (GQ477367); FIPV, feline infectious 
peritonitis virus (AY994055); MCoV, mink CoV (HM245925); RhBatCoV HKU2, rhinolophus bat CoV HKU2 (EF203064); Mi-BatCoV 1A, 
miniopterus bat CoV 1A (NC_010437); Mi-BatCoV 1B, miniopterus bat CoV 1B (NC_010436); Mi-BatCoV HKU8, miniopterus bat CoV 
HKU8 (NC_010438); HCoV-229E, human CoV 229E (NC_002645); PEDV, porcine epidemic diarrhea virus (NC_003436); Sc-BatCoV-512, 
scotophilus bat CoV 512 (NC_009657); HCoV-NL63, human CoV NL63 (NC_005831); Hi-BatCoV HKU10, hipposideros bat CoV HKU10 
(JQ989266); Ro-BatCoV HKU10, rousettus bat CoV HKU10 (JQ989270); RbCoV HKU14, rabbit CoV HKU14 (JN874559); HCoV-OC43, 
human CoV OC43 (NC_005147); PHEV, porcine hemagglutinating encephalomyelitis virus (NC_007732); CRCoV, canine respiratory 
CoV (JX860640); BCoV, bovine CoV (NC_003045); GiCoV, giraffe CoV (EF424622); SACoV, sable antelope CoV (EF424621); ECoV, 
equine CoV (NC_010327); HCoV-HKU1, human CoV HKU1 (NC_006577); MHV, murine hepatitis virus (NC_001846); RCoV, rat CoV 
(NC_012936); SARSr-CiCoV, SARS-related palm civet CoV (AY304488); SARS-CoV, severe acute respiratory syndrome—associated 
human CoV (NC_004718); SARSrCoV CFB, SARS-related Chinese ferret badger CoV (AY545919); SARSr-Rh-BatCoV HKU3, SARS- 
related rhinolophus bat CoV HKU3 (DQ022305); MERS-CoV, Middle East respiratory syndrome CoV (JX869059); Pi-BatCoV HKU5, 
pipistrellus bat CoV HKU5 (NC_009020); Ty-BatCoV HKU4, tylonycteris bat CoV HKU4 (NC_009019); Ro-BatCoV HKU9, rousettus bat 
CoV HKU9 (NC_009021); IBV, infectious bronchitis virus (NC_001451); IBV-partridge, partridge CoV (AY646283); TCoV, turkey CoV 
(NC_010800); IBV-peafowl, peafowl CoV (AY641576); BWCoV-SW1, Beluga whale CoV SW1 (NC_010646); BdCoV HKU22, bottlenose 
dolphin CoV HKU22 (KF793824); WiCoV HKU20, wigeon CoV HKU20 (JQ065048); CMCoV HKU21, common-moorhen CoV HKU21 
(JQ065049); NHCoV HKU19, night-heron CoV HKU19 (JQ065047); BuCoV HKU11, bulbul CoV HKU11 (FJ376619); MunCoV HKU13, 
munia CoV HKU13 (FJ376622); MRCoV HKU18, magpie—robin CoV HKU18 (JQ065046); ThCoV HKU 12, thrush CoV HKU12 (FJ376621); 
WECoV HKU16, white-eye CoV HKU16 (JQ065044); PorCoV HKU15, porcine CoV HKU15 (JQ065042); SpCoV HKU17, sparrow CoV 
HKU17 (JQ065045). Numbers at nodes represent bootstrap values. Scale bar indicates the estimated number of substitutions per 5 aa. 
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Figure 6. Sodium dodecyl sulfate—polyacrylamide gel electrophoresis 
and Western blot analysis of a novel coronavirus, dromedary camel 
coronavirus UAE-HKU23, discovered in dromedaries of the Middle 
East, 2013. Nucleocapsid protein was expressed in Escherichia coli. 
M, protein molecular-mass marker; kDa, kilodaltons. Lanes: 1, non- 
induced crude E. coli cell lysate; 2, induced crude E. coli cell lysate 
of DcCoV UAE-HKU23 nucleocapsid protein; 3, purified recombinant 
DcCoV UAE-HKU23 nucleocapsid protein; 4, dromedary camel 
serum sample strongly positive for antibody against nucleocapsid 
protein; 5, dromedary camel serum sample moderately positive for 
antibody against nucleocapsid protein; 6 and 7: dromedary camel 
serum sample negative for antibody against nucleocapsid protein. 


The Ka/Ks of all the coding regions in DcCoV UAE- 
HKU23 was <0.5. 

By using the uncorrelated relaxed clock model on 
RdRp gene sequences, we estimated the date of diver- 
gence between DcCoV UAE-HKU23 and BCoV to be ~46 
years ago. We estimated that the 3 strains of DcCoV UAE- 
HKU23 diverged from their most recent common ancestor 
in March 2010 (the 95% highest posterior density interval, 
August 2006—September 2012) (online Technical Appen- 
dix Figure 1). 


Discussion 

We discovered a novel CoV, but no MERS-CoV, in 
dromedaries from the Middle East. Dromedaries are | of 
2 surviving camel species. Dromedaries (C. dromedarius; 
1-humped camels) inhabit the Middle East and northern 
and northeastern Africa; Bactrian camels (C. bactrianus, 
2-humped camels) inhabit Central Asia. Among the 20 mil- 
lion camels on earth, 90% are dromedaries. In 2012, there 
were ~360,000 dromedaries in the United Arab Emirates. 

In this study, we discovered a novel CoV, DcCoV 
UAE-HKU23, from 4.8% of 293 fecal samples collected 
from dromedaries in Dubai. The positive samples were not 
collected from the same farm or stable. Moreover, there 
was >0.2% nt difference among the 3 complete genomes 
sequenced, indicating that the positive samples were not 
collected from a clonal outbreak. In our study, 21.3% of 
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dromedary calves, but only 0.4% of adult dromedaries, 
were RT-PCR positive for DcCoV UAE-HKU23; this 
finding indicates that dromedary calves are probably more 
susceptible than adult dromedaries to infection with Dc- 
CoV UAE-HKU23. Furthermore, DcCoV UAE-HKU23 is 
probably stably evolving in dromedaries because the Ka/ 
Ks of all the coding regions in the genome were <0.5. In 
this study, 4 of the 12 positive samples were collected from 
dromedaries with diarrhea. A previous report also described 
the presence of a betacoronavirus in the fecal sample of a 
dromedary calf with diarrhea (35). This finding raises the 
question of the pathologic significance of DcCoV UAE- 
HKU23 for camelids and warrants further animal studies. 

Our serologic data showed little cross-reactivity be- 
tween DcCoV UAE-HKU23 and SARS-CoV, Pi-BatCoV 
HKUS, and Ro-BatCoV HKU9. This finding is in line 
with findings from our previous studies of Ro-BatCoV 
HKU9, which also showed minimal serologic cross-re- 
activity among the 4 lineages of betacoronaviruses (/6). 
These results suggest that there should be minimal cross- 
reactivity between DcCoV UAE-HKU23 and MERS- 
CoV, which belong to 2 different CoV lineages. Because 
we showed an extremely high prevalence of MERS-CoV 
antibodies in the serum samples by Western blot analy- 
sis, indirect immunofluorescence, and neutralization anti- 
body testing, concurring with findings in a previous study 
(24), we would also expect a similar high prevalence of 
DcCoV UAE-HKU23 antibodies if there was major se- 
rologic cross-reactivity between MERS-CoV and DcCoV 
UAE-HKU23. However, our serologic data only revealed 
the presence of DcCoV UAE-HKU23 antibodies in 52% 
of the serum samples, indicating that no correlation exists 
between seropositivity to DcCoV UAE-HKU23 and sero- 
positivity to MERS-CoV. Furthermore, we found no cor- 
relation between seropositivity to DcCoV UAE-HKU23 
and MERS-CoV antibody titers. 

In this study, correlation between DcCoV UAE- 
HKU23 RT-PCR positivity and seropositivity also 


Table 4. Detection of antibodies to MERS-CoV in dromedaries in 
the Middle East, 2013* 

Test, antibody titer 

Indirect immunofluorescence 


No. (%) positive samples 


160 2 (3.4) 
320 12 (20.3) 
640 16 (27.1) 
1280 11 (18.6) 
2560 14 (23.7) 
5120 4 (6.8) 
Neutralization antibody test 
<10 1 (1.7) 
10 7 (11.9) 
20 14 (23.7) 
40 28 (47.5) 
80 8 (13.6) 
160 1 (1.7) 


*MERS-CoV, Middle East respiratory syndrome coronavirus. 
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Table 5. Estimates of nonsynonymous and synonymous 
substitution rates in the genomes of a novel betacoronavirus, 
DcCoV UAE-HKU23, discovered in dromedaries of the Middle 
East, 2013* 


Gene Ka Ks KalKs 
NSP1 0 0.004 0 
NSP2 0 0.002 0 
NSP3 0.001 0.006 0.167 
NSP4 0.001 0.007 0.143 
NSP5 0.001 0.003 0.333 
NSP6 0.001 0.003 0.333 
NSP7 0 0 - 
NSP8 0.016 0 
NSP9 0) 0 - 
NSP10 0 0) - 
NSP11 0 0 - 
NSP12 0 0.002 0 
NSP13 0 0.002 0 
NSP14 0.001 0 - 
NSP15 0.001 0 - 
NSP16 0 0 - 
NS2 0.001 0.008 0.125 
HE 0.001 0.002 0.5 
Spike 0.001 0.004 0.25 
NS5 0 0.009 0 
Envelope 0 0 - 
Membrane 0 0 - 
Nucleocapsid 0 0 - 
N2 0 0 - 


DcCoV, dromedary camel coronavirus; Ka, nonsynonymous site; Ks, 
synonymous site; NSP, nonstructural protein; NS, nonstructural; HE, 
hemagglutinin esterase; N, nucleocapsid. 


cannot be ascertained because the fecal samples and se- 
rum samples were collected from different dromedaries. 
Because MERS-CoV was not present in dromedaries 
in the present study, an intensive search in dromedar- 
ies and other animals in other locations in the Middle 
East would be helpful in the search for the animal source 
of MERS-CoV. 


—, 
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DcCoV UAE-HKU23 is a member of betacoronavirus 
Al (Figure 7). Comparison of the amino acid identities of 
the 7 conserved replicase domains for species demarcation 
(i.e., ADP-ribose 1”-phosphatase, NSP5 [3CL?°], NSP12 
[RdRp], NSP13 [helicase], NSP14 [ExoN], NSP15 [Ne- 
ndoU], and NSP16 [O-MT]) (36) between DcCoV UAE- 
HKU23 and other CoVs of betacoronavirus Al revealed 
that in all 7 domains, the amino acid sequences of DcCoV 
UAE-HKU23 and other betacoronavirus Al members 
shared >90% identity. This finding indicates that DcCoV 
UAE-HKU23 should be a member of betacoronavirus Al. 

Furthermore, the genome characteristics of DcCoV 
UAE-HKU23 showed features similar to those of other 
betacoronavirus Al members. The genomes of all beta- 
coronavirus Al members have G+C contents of 0.37 and 
the genome of RbCoV HKU 1/4, a recently discovered CoV 
closely related to betacoronavirus Al (/9), has a G+C con- 
tent of 0.38. This G+C content differs substantially from 
those of other CoV species of lineage A betacoronavi- 
ruses, which have G+C contents of 0.32 (HCoV-HKU1), 
0.41 (rat CoV) and 0.42 (murine hepatitis virus) (Table 2). 
The difference in the genome characteristics between Al 
and non-Al members of betacoronavirus is also reflected 
by their codon usage bias, in which the general preference 
of using G/C in the third position of the codons decreases 
from murine hepatitis virus and rat CoV to betacoronavirus 
Al members and RbCoV HKU14 to HCoV-HKUI1 (online 
Technical Appendix Figure 2). The cleavage site for spike 
protein of betacoronavirus Al members is RRS/QRR, 
whereas those of HCoV-HKUI, RbCoV HKU14, and 
murine hepatitis virus are RRKRR, LRSRR, and RAR/H/ 
DR/S, respectively. The length of membrane genes for 
betacoronavirus Al members and RbCoV HKU14 is 693 
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bases, whereas the lengths for HCoV-HKU1 and murine 
hepatitis virus are 672 and 687 bases, respectively. 

DcCoV UAE-HKU23 is phylogenetically closely re- 
lated to other betacoronavirus Al members, which in turn 
are closely related to other CoVs of betacoronavirus lin- 
eage A. Despite their close relationships, no recombination 
was detected between DcCoV UAE-HKU23 and other be- 
tacoronavirus Al members by bootscan analysis (data not 
shown). These CoVs of betacoronavirus Al may be using 
different receptors in their corresponding hosts because 
their spike protein is one of the proteins that show the larg- 
est difference among different CoVs. Most of the differ- 
ences among the spike proteins in different betacoronavi- 
rus Al members were also observed in the N terminal half 
of their spike protein, where the receptor binding domains 
should be located. 

Camels are one of the most unique mammals on 
earth. In particular, they have shown perfect adaptation 
to desert life, which presents temperature extremes and 
a scarce supply of food and water. In the past, camels 
were used for transportation of humans and goods and for 
military uses. Moreover, for humans, they provide a good 
source of meat, milk, and wool. Camels are also impor- 
tant recreational animals in the Middle East and are used 
for camel racing. Having been associated with humans for 
at least 5,000 years, camels usually pose little physical 
danger to humans. However, infectious pathogens, such 
as brucellosis, can occasionally be transmitted from cam- 
els to humans. Apart from the present novel CoV, viruses 
of at least 8 taxonomic families (i.e., Paramyxoviridae, 
Flaviviridae, Herpesviridae, Papillomaviridae, Picorna- 
viridae, Poxviridae, Reoviridae, and Rhabdoviridae) have 
been found to infect camels (37-39). Because camels are 
closely associated with humans, continuous surveillance 
of viruses in this hardy group of animals is needed to un- 
derstand the potential for virus emergence and transmis- 
sion to humans. 
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